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Abstract. As a technical supplement to the above mentioned paper we 
present 192 consecutive eigenstates for the Robnik billiard with the shape 
parameter A = 0.15 from 10,001st to 10,192nd, by showing the plots in the 
configuration space and in the phase space. The latter is smoothed pro- 
jection of the Wigner function onto the surface of section. By comparison 
with the classical SOS plots we thus examine all eigenstates and classify 
them in regular and irregular: There are 70 regular states and 122 irregular 
states, thus giving the estimate of the relative measure of the regular com- 
ponent pi = 0.365, which is in excellent agreement with the classical value 
pi = 0.360 calculated and reported by Prosen and Robnik (1993). 
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Explanation of plots 



In this work we provide additional material to our recent paper (Li and Rob- 

nik 1994a). The method and the mathematics of analyzing the eigenstates 
in the Robnik billiard defined as the quadratic complex conformal map of 
the unit disk (in the z-plane, \z\ < 1) onto the complex w-plane, namely 

Bx^ {w\w ^ z + \z^, 1^1 < 1}. (1) 
with A = 0.15, is presented in our recent paper (Li and Robnik 1994b). 

The wavefunctions we are looking at are the eigenfunctions of the Schrodinger 
equation (Helmholtz equation): 

+ = 0, * = at the boundary of Bx (2) 

where E = k"^ is the eigenenergy and k the wavenumber. (So we are using 
units such that Planck's constant h = 1 and 2m = 1, where m is the mass of 
the point billiard particle.) 

In the configurational plots in figures 1-8 we show the contours of con- 
stant probability density for the eigenstates 10,001st through 10,192nd with 
each figure containing 24 consecutive plots ordered as left-right top-down se- 
quence. The scale and the coordinates are uniquely specified by the equation 
defining the boundary namely, w — z + Xz'^, where \z\ = 1. The contours of 
constant probability density are chosen in steps of 1/8 of the maximal value 
of each individual plot. 

In order to investigate the eigenstates in the quantum (Wigner) phase space 
we have first to define the classical phase space and the surface of section. 
The usual bounce map (Poincare map) in the Birkhoff coordinates (arclength 
versus tangent unit velocity vector component) is not suitable for our pur- 
pose, because \t' vanishes on the boundary. Therefore we choose the surface 
of section defined hj v = lm{w) = 0: Our surface of section is now specified 
by the crossing point coordinate u on the abscissa versus the conjugate mo- 
mentum equal to the tangential component of the velocity vector of length 
k with respect to the line of section v = 0. In figure 9 we show 24 identical 
plots of the geometry of the largest chaotic component for A = 0.15, in ex- 
actly the same size as the quantal phase space plots which we will show in 
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figures 10-17, which enables the reader to perform the comparison of classical 
and quantal plots and the classification of the eigenstates in regular and ir- 
regular, by overlaying the sheets. We do not show further details of the KAM 
scenario inside the stability islands in order not to obscure the structure of 
the phase space. 

The Wigner function (of an eigenstate f)) defined in the full phase 
space {u,v,pu,Pv) is 

W{q, p) = ^ J d'y. exp(-2p ■ X)vl/(q - X/2)vl/(q + X/2) (3) 

where we have specialized to our real case, and also two degrees of freedom 
and ^ = 1. Here q = {u,v) and p = {pu,Pv)- In order to compare the 
quantum Wigner functions with the classical Poincare maps on the surface 
of section we define the following projection of (H) given as 

psos{u,Pu) = J dpyW{u,0,pu,Pv), (4) 
which nicely reduces the number of integrations by one and is equal to 

1 /" XX 

Psos{u,Pu) = c?a;exp(ia;p„)^(tt + -,0)^(m - -,0) (5) 
/vr J z z 

As is well known the Wigner function and its projections are not positive def- 
inite and indeed one typically finds small and inconvenient but nevertheless 
physical oscillations around zero which seriously obscure the main structural 
features. Therefore in order to compare the classical and quantal phase space 
structure it is advisable to smooth the Wigner function or its projections (|^) 
by a normalized Gaussian kernel with a suitably adapted dispersion. Such 
procedure has been introduced and used in (Takahashi 1989, Leboeuf and 
Saraceno 1990, Heller 1991, Prosen and Robnik 1993), which is Husimi type 
representation but the effective area of our Gaussian kernel will be smaller 
than 2n. 

In figures 10-17 we show the smoothed object for each consecutive eigen- 
state from the 10,001st to the 10,192nd, each figure containing 24 states 
arranged in order left-right top-down. The lowest contour is at the level of 
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0.15 of the maximal value and the step size upwards is also 0.15 of the max- 
imum. 

The examination of the states in figures 10-17 and classification leads to the 
following result: There are 70 regular and 122 irregular states. This yields 
the relative fraction of regular levels pi = 0.365, in excellent agreement with 
the classical value pi = 0.360 (Proscn and Robnik 1993). 
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